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2.4 Spreading of a shallow mass on an incline 
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The spreading of a finite mass of paint, paper pulp, mud or lava on an inclined plane is of 
interest to a variety of industrial and geological problems. For mud modeled as a Bingham- 
plastic non-Newtonian fluid, Liu & Mei (1989) solved the equation similar to (2.3.22) numer- 
ically. An analytical solution for Herschel-Bulkley fluid was given later by Huang & Garcia 
(1997). We modify their theories for the simpler case of Newtonion fluid here. 

The free surface of a thin layer is expected to flatten in time over most of the profile, 
but it should steepen near the downstream front where the spatial derivative is much more 
important than elsewhere. Let us study the problem by dividing the total fluid extent into 
two: the far field not too close to the steep front and the near field around the front. 

2.4.1 Far field away from the front 

Let the total length of the thin layer be L and the maximum depth be H, with H/L <^ 1. 
We introduce the normalized variables suitable for the far field. 



h = Hh', X = Lx', t = Tt' 



(2.4.1) 



where T will be chosen to simplify the appearance of the final equation. Thus, 
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Let us choose the time time scale to be 
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(2.4.2) 
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where the denominator is the average fluid speed across the depth. The dimensionless 
equation for the farfield reads, 



dh!_ d_ 
dt' dx' 
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= (2.4.3) 



Assume 

^ = e < 1 (2.4.4) 
Ltan^ ^ ^ 

(2.4.3) can then be approximated by the hyperbohc equation, 

dh' dh''^ dh' , ,odh' ^ , , , 

It is of the class called kinematic wave equation in flood hydrology. Solution can be obtained 
by the theory of characteristics. 

In particular (2.4.5) can be rewritten in the form 

Oh Oh 

—dt + —dx = d/i = (2.4.6) 
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if 

^ = 3h'^ (2.4.7) 

dt' ^ ' 

The last two equations imply that h' remains constant for all t' along the characteristic 

curve ,T'(t') defined by the differential equation (2.4.7). Moreover, if the initial profile is 
presecribed, 

U{x\^) = UXx), (2.4.8) 

then all characteristics are straight but of different slopes. The characteristic originated from 
the initial point x' = ^' at t' = is the straight line 

^(^,0 = 3/i^(0i + e, (2.4.9) 

along which 

/iv,o = ^;(o (2.4.10) 

In principle we can solve for ^ from (2.4.9) in term of and substitute the result 
into (2.4.10) to get h'{x',t'). For any initial hump, the characteristics at the front must 
intersect one another, impling multivaluedness of solution, which is phycially unacceptable. 
The solution can still be proceeded if a discontinuity, shock, is allowed at x' — x'g{t'), as long 
as mass is conserved. 

As an example, consider a triangular initial profile: 
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The profile has a shock front to begin with, which is a mathematical idealization of course. 
Prom (2.4.9) we get 

1 c', r<o,r>i 
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Within < ^' < 1, ^' can be solved in terms of x',t', 

r = ^ (-1 + 71+12^) 

For any t' > 0, 



and 



for x' <0, k x' > x'^{t'), 
/i' = e' = ^(-l + Vl + 12a;'t'), 0<x'<x'^{t). 



The shock front x'g{t') is unknown. 

To locate the shock front x'g{t'), we invoke mass conservation, 

1 f-^'si^') i< dx' 
initial volume — - — h' dx' = h' — - 
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Prom (2.4.12), 
so that 

Thus 
or 



x' = 3h'H' + h' 
= 6h't' + 1 
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(2.4.12) 



(2.4.13) 

(2.4.14) 
(2.4.15) 



(2.4.16) 



(2.4.17) 



(2.4.18) 



This gives the shock height h'^ implicitly as a function of t' . To get the shock postion x'g{t') 
we put this result into (2.4.17) 
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(2.4.19) 



Together (2.4.18) and (2.4.19) give x'g{t'), as plotted in figure (2.4.1). To be needed later we 
note that the shock speed is given by 



dx'g _ dx'g dt' _ 1 - ?,h'~'^ 4/i'/ 
ItF ~ lih/'dh' ~ 4 1 - 3/i': 
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(2.4.20) 



Figure 2.4.1: Time variation of the location x'g{t') and depth h'^{t') of shock front. 



2.4.2 Near field of downstream front 



Clearly near the shock the local free surface slope cannot satisfy (2.4.4). We define the 
longitudinal length scale for the near field to be Lg so that 



(2.4.21) 



and renormalize the near field, 



h = HH', t = Tt', ^^M = x', or x' = x', + eX' (2.4.22) 

Thus the near field solution H' is a function of X'{x',t') and t'. The derivatives are now 
related by the chain rule, 



dH'{x', t') _ dH' dH' dX' _dH' _1 dx'^ dH' 
dt' ~ 'df ^ dX'~dF' ~ 'W ^ 7 dt' dX' 

dH'{x',t') dH' dX' IdH' 



Eq. (2.4.3) becomes 
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To the leading order the near field is approximately governed by 
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(2.4.23) 
(2.4.24) 

(2.4.25) 
(2.4.26) 
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Figure 2.4.2: Matching the near-field and far-field at the front by preserving mass. 



which is the ordinary differential equation for a stationary wave, moving at the known shock 
speed dx'g/dt' = h'^ . Integrating with respect to X' once, we get 



0. 



The constant of integration is taken to be zero so that the surface touches the dry bed where 
H' = 0. We may rewrite (2.4.26) as 



dX' 



H'^dH' 



dH' 
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It follows upon integration that 



X'-X' 



(2.4.27) 



This is a smooth surface decreasing monotonically from H' = h'^ at X' = — oo to H' = 
(dry bed) at X' = X'^. To determine the value of X^:, we require that mass under the smooth 
profile be the same as that under the profile with a shock, as sketched in Figure 2.4.2. 



{h'^ - H') dX' 



XL 



H\X')dX' 



Referring to Figure ??, the integrals can be equivalently carried out as 



H'(0) 



X\H')dH' 



X\H')dH' 



where H'{0) is the height at X' = 0, given by 

fK-H'joy 

H'{0), 
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(2.4.28) 



(2.4.29) 



(2.4.30) 
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Substituting (2.4.27) into (2.4.29), we get, for the left-hand-side. 



and for the right-hand side. 



Therefore, 



K f h' - H'\ 
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1 /l-^ 

log 7^ 
\1 + z 



dz + 1 



Using the facts 



we get 



log(l - z)dz = -1, log(l + z)dz = 2 log2 - 1 
Jo Jo 

Q-log2) = 0.19897/^; 

The near field solution is now complete. 

Recall that the near and far fields coordinates are related by 



(2.4.31) 



X = x„ + eX 



in general and 



<(0=^'.(0 + ^Q-log2) Kit') 



(2.4.32) 



(2.4.33) 



in particular. 

We can get the uniformly valid solution (hu.v.) by adding the near and far field solutions 
and subtracting the common part. Thus, in terms of the far-field variables. 



h'u,v. = h' + H'- K 



(2.4.34) 



where h' is given by (2.4.15) and H' is given by (2.4.27), while the coordinates are related 
by (2.4.32). 

Sample profiles of the free surface are shown for e = 0.1 in Figure 2.4.3 for various 
instants. 
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Figure 2.4.3: Slow evolution of the free surface of a vicous fluid down a dry incline. The 
initial profile is triangular with e = 0.1. 



